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Preface 

This  report  is  our  attempt  to  determine  the  feasibility 
of  an  earth-moon  mission  using  a  recursive  navigation  technique 
where  only  individual  earth-moon  fixes  are  used  for  navigational 
measurements.  We  believe  that  the  report  is  organized  in  such  a 
manner  as  to  permit  easy  assimulation  of  the  theories  and  techniques 
by  the  reader  with  some  experience  in  lunar  navigation  studies. 
However,  we  do  hope  that  the  uninitiated  can  glean  some  informative 
meaning  of  the  results  and  conclusions  recorded.  Since  the 
introduction  is  the  key  to  our  pattern  of  development  and  sets  forth 
the  assumptions  upon  which  this  study  is  based,  the  reader  should 
digest  it  before  attempting  the  body  of  the  report. 

This  study  is  by  no  means  independent  in  the  true  sense  of  the 
word.  As  is  obvious  from  the  bibliography,  it  is  based  on  the 
findings  of  several  experts  in  the  field  of  astronautical  guidance. 
In  particular,  we  would  like  to  single  out  the  work  of  Dr.  Richard. 

H.  Battin  of  the  M.I.T.  Instrumentation  Laboratory.  Dr.  Battin’s 
well -formulated  recursive  navigation  theory  and  matrix  techniques 
are  used  exclusively  in  this  study.  It  is  safe  to  say  that  without 
them  this  study  would  have  never  been  attempted.  Also,  we  must 
acknowledge  our  faculty  thesis  advisor,  Capt.  Charles  J.  Januska, 
for  suggesting  this  study  and  giving  us  the  benefit  of  his  insight 
into  Dr.  Battin* s  theory.  Capt  Januska* s  constant  advice  and 
encouragement  made  it  possible  for  us  to  complete  this  study.  Last, 
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but  by  no  means  least,  we  are  grateful  for  the  I.B.M.  7090 
digital  computer  programming  assistance  of  Mr,  William  R.  Poole, 
a  mathematician  in  the  A.S.D.  Digital  Computation  Division  (WPAFB, 
Ohio).  Without  Mr.  Poole’s  capable  programming  of  our  navigation 
equations  and  his  interest  and  patience  in  making  the  many  changes 
suggested  as  the  study  progressed,  we  could  have  never  completed 
this  study  on  time. 

This  investigation  has  been  of  great  benefit  to  us  in  two 
distinct  ways.  First,  we  have  become  familiar  with  space  navigation 
theories  and  gained  an  insight  into  the  statistical  theory  upon 
which  they  are  based.  Secondly,  we  have  gained  invaluable  experience 
in  digital  programming  techniques. 


Paul  J.  Ellmer 
Donald  R.  Hefty 
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Abstract 

The  purpose  of  this  study  is  to  determine  the  feasibility 
of  using  a  recursive  navigation  theory  in  which  individual  earth- 
moon  angular  measurements  (ffoces)  are  processed  as  they  are  made 
and  combined  with  the  current  best  estimate  of  position  and  velocity 
to  produce  an  improved  estimate.  The  recursive  navigation  techniques 
of  Dr.  Battin  of  the  M.I.T.  Instrumentation  Laboratory  are  used. 
Estimates  made  with  an  earth-moon  fix  interval  of  5  minutes  or 
less  are  compared  to  estimates  made  with  conventional  earth-star, 
moon-star  fixes  of  15  minutes.  Results  show  the  earth-moon  fix 
to  be  a  feasible  technique. 
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UNCERTAINITY  IN  THE  ESTIMATE  OF  POSITION  AND  VELOCITY 
ALONG  A  LUNAR  TRAJECTORY  USING  INDIVIDUAL 
EARTH-MOON  ANGULAR  MEASUREMENTS  IN  A 
RECURSIVE  NAVIGATION  THEORY 


I.  Introduction 


The  design  and  development  of  self-contained  lunar  space 
navigation  techniques  are  of  paramount  interest  to  many  of  the 
space  scientists  and  astronautical  engineers  of  today.  The  most 
critical  factors  in  the  design  of  such  techniques  are  the  accuracy 
and  reliability  of  the  spacecraft  equipment  needed  to  implement 
them.  Since  the  navigation  system  is  one  of  the  most  complex  of 
the  spacecraft  subsystems,  special  attention  has  been  given  to 
reliability.  Generally,  this  reliability  is  enhanced  by  the  use 
of  simple  equipment  concepts.  Therefore,  efforts  to  design  new 
techniques,  as  well  as  modify  present  techniques,  utilizing  simpler 
equipment  concepts  is  a  continuous  process.  This  study  is  concerned 
with  the  simplification  of  a  proven  concept  of  space  navigation  as 
applied  to  an  earth-moon  mission  and  is  justified  on  the  basis  that 
a  simplification  of  equipment  can  be  obtained  without  sacrifice  to 
accuracy. 

Statement  and  Analysis  of  the  Problem 

A  unique  injection  ( launch) velocity  imparted  to  an  earth  orbiting 
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vehicle  at  a  specific  point  in  space  will  result  in  an  unique 
preplanned  trajectory.  As  s  ~esult  of  inject .on  velocity  errors, 
however,  the  spacecraft  will  not  follow  this  preplanned  trajectory. 
Therefore,  in  order  to  determine  the  actual  trajectory,  navigational 
information  must  be  obtained  during  the  flight.  A  source  of  this 
navigational  information  in  space  is  the  observation  of  the  angle 
between  selected  celestial  bodies  as  seen  from  the  spacecraft.  It 
can  be  readily  appreciated  that  an  angular  measurement  between  two 
celestial  bodies  is  unique  to  only  one  point  in  space  for  any  given 
time.  Each  angular  measurement  (navigational  measurement,  naviga¬ 
tional  fix,  or  simply,  fix)  need  not  determine  the  position  and 
velocity  of  the  spacecraft;  the  information  from  a  series  of 
individual  fixes  taken  over  a  period  of  time  may  be  combined  as  a 
recursive  operation  in  which  individual  measuruments  are  processed 
as  they  are  made  and  combined  with  the  current  beat  estimate  to 
produce  an  improved  estimate.  Theoretically,  six  such  fixes  would 
suffice  to  determine  the  position  and  velocity  of  the  spacecraft. 
However,  due  to  the  inherent  inaccuracies  of  the  onboard  angular 
measurement  devices  (an  optical  sextant  in  this  study),  statistical 
methods  baaed  on  the  known  statistical  behavior  of  these  devices 
are  necessary  to  ascertain  the  uncertainty  (r.m.s.  error)  in  the 
position  and  velocity  estimate  obtained. 

The  recursive  navigation  the  ary  of  Dr.  Richard  H.  Battir  of 
the  M.I.T*  Instrumentation  laboratory  is  used  in  this  study 
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(Ref  3:303-340).  This  theory  involves  the  formulation  of  an 
optimum  linear  estimate  as  a  recursive  operation  in  which  the 
best  estimate  of  position  and  velocity  obtained  by  the  beat  choice 
of  navigational  fix  is  combined  with  newly  acquired  information  to 
produce  an  improved  estimate.  The  derivation  of  the  optimum  linear 
operation  is  baaed  on  the  application  of  a  simple  leaat-squares 
estimation  technique. 

The  Problem.  This  study  is  concerned  with  the  determination 
of  the  uncertainties  in  the  estimates  of  position  and  velocity 
along  a  lunar  trajectory  using  individual  earth -moon  angular 
measurements  in  a  recurjdve  navigation  theory.  The  -unar  trajectory 
in  this  study  is  define!  as  an  elliptical  path  from  a  parking 
orbit  about  the  earth  to  a  point  near  the  moon's  sphere  of  influence 
(see  Appendix  A). 

In  previous  studi.es  made  with  the  recursive  navigation  theory 
outlined,  navigational  measurements  were  made  along  the  spacecraft's 
trajectory  using  various  star-planet  combinations.  At  each  pre¬ 
planned  measurement  time,  the  particular  star  and  planet  (in  a  lunar 
trajectory  study,  several  prominent  stars  and  the  planets,  earth 
and  moons  are  usually  considered)  contributing  the  least  uncertainty 
in  the  estimate  of  position  aDd  velocity  is  selected.  The  analysis 
of  the  studies  using  this  method,  referred  to  as  Fix  2  in  this  study, 
shows  that  acceptable  accuracy  was  obtained  (Ref  2:35~3^)» 

In  this  study  the  fix  at  each  preplanned  point  in  time  ilong 
the  spacecraft's  trajectory  was  restricted  to  the  measurement  of 
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the  angle  between  the  earth  and  moon  only.  This  method  is  referred 
to  as  Fix  1  in  this  study.  At  first  glance,  this  restriction  seems 
to  be  unworthy  of  investigation  because,  given  the  same  fix  times 
for  Fix  1  and  Fix  2,  it  is  obvious  that  Fix  2  is  superior  to  Fix  1. 
However,  considering  smaller  time  intervals  between  the  measure¬ 
ments  of  Fix  1,  the  situation  takes  on  a  new  prospective.  In  Dr. 
Battin's  studies.  Fix  2  (Ref  2:28)  was  restricted  to  time  intervals  of 
no  l*»ss  than  15  minutes  because  of  the  time  required  to  reposition 
the  optical  sextant  to  the  next  star-planet  fix.  In  Fix  1,  since 
the  sextant  is  measuring  only  the  earth-moon  angle,  no  complex 
positioning  arrangement  need  be  implemented.  Since  the  earth  and 
moon  is  continuously  monitored,  fixes  can  be  taken  at  much  smaller 
time  intervals.  This  allows  more  measurements  to  be  made  over  a 
given  span  of  time  and  increases  the  accuracy  of  the  estimate  of 
position  and  velocity  at  the  end  of  that  span . 

Time  intervals  as  small  as  30  seconds  were  considered  practical 
for  Fix  1,  since  the  only  tine  restriction  of  any  consequence  was 
considered  the  processing  time  of  the  spacecraft’s  onboard  digital 
computer.  Therefore,  in  essence,  this  study  is  an  investigation 
of  using  Fix  1  along  a  lunar  trajectory  at  time  intervals  less 
than  15  minutes  in  a  recursive  navigation  theory  to  determine  whether 
or  not  it  will  give  results  comparable  to  Fix  2  at  a  time  interval 
of  15  minutes. 

Pattern  of  Development .  The  following  sections  attempt  to  outline 
the  recursive  navigation  theory  and  show  how  it  is  applied  to  Fix  1 
and  Fix  2  to  obtain  the  uncertainties  in  the  estimate  of  position 
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and  Telocity  along  a  lunar  trajectory. 

Section  II  contains  a  brief  outline  of  the  nature  of  the 
single  observation  fix  and  the  linearized  trajectory  equations 
which  are  used  in  the  recursive  navigation  theory.  This  section 
is  supplemented  by  Appendixes  A,  B,  and  C.  Appendix  A  describes 
the  basic  considerations  in  planning  a  lunar  trajectory.  Bie 
navigational  fix  equations  are  presented  in  Appendix  B,  while 
Appendix  C  presents  the  parameters  of  the  lunar  trajectory  used 
in  this  study  and  states  the  equations  which  represent  the 
deviations  in  position  and  velocity  along  the  reference  trajectory. 

Section  III  gives  some  basic  statistical  considerations 
for  the  navigation  theory  and  outlines  the  equations  which  define 
the  uncertainties  in  the  estimate  of  position  and  velocity  at 
any  point  in  tims  along  the  spacecraft's  trajectory.  Appendix 
D  supplements  Section  III  by  presenting  the  assumed  initial 
errors  at  injection  of  the  spacecraft  into  its  lunar  trajectory 
and  the  assumed  error  of  the  optical  sextant  considered. 

The  results  of  Section  IY  graphically  portray  the  comparison 
between  Fix  1  ana  Fix  2  for  different  time  intervals  along  the 
trajectory.  In  addition  to  this  comparison »  other  factors  observed 
which  relate  to  the  uncertainty  in  the  position  and  velocity 
estimate  using  - • x  1  are  presented.  These  factors  are: 

1.  The  effect  of  the  inclination  between  the  plane 
of  the  spacecraft's  reference  (preplanned)  trajectory  and  the 
orbital  plane  of  the  moon. 
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2.  The  effect  of  the  assumed  r.m.s.  error  of  the 
onboard  optical  sextant . 

3.  The  effect  of  including  in  Fix  l  the  apparent  angular 
diameter  measurements  of  the  earth  and  moon. 

The  conclusions  of  Section  V  discuss  the  results  obtained  using 
Fix  1  and  their  meaning  in  terms  of  the  compatability  with  uncertain¬ 
ties  in  position  and  velocity  using  Fix  2.  Possible  extensions  to 
this  study  are  outlined  in  the  recommendations  of  Section  VI. 

Assumptions .  The  solution  of  this  problem  was  based  on  the 
following  assumptions: 

1.  The  spacecraft's  motion  is  adequately  described  in  a  single 
inverse  squared  gravity  field  (earth).  This  assumption  is  discussed 
in  more  detail  in  Appendix  A. 

2.  The  location  of  stars  was  arbitrarily  selected  along 
the  six  axes  of  the  space  coordinate  system  (Appendix  C)  and  at 

an  infinite  distance.  This  assumption  is  valid  because  the 
direction  cf  the  stars,  corrected  for  aberration  and  parallax  for 
the  period  of  a  lunar  mission,  are,  for  all  practical  purposes, 
unvarying  regardless  of  vehicle  position.  Also,  aboard  the  space¬ 
craft  the  stars  are  unobscured  and  appear  as  point  sources  of  light. 

3.  The  features  of  the  moon  can  be  readily  observed 
from  a  spacecraft  and  the  terrain  features  on  earth  will  also  be 
observable  from  the  spacecraft  despite  some  degradation  due  to 
atmospheric  effects.  Therefore,  the  use  of  an  onboard  optical 
sextant  suggests  itself  as  one  possible  technique  for  making  angular 
measurements . 
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4.  The  only  equipment  inaccuracies  considered  are  in 
the  optical  sextant  of  the  spacecraft’s  navigational  system.  In 
addition,  the  optical  sextant  is  assumed  to  have  a  known  random 
and  independent  error.  Section  III  and  Appendix  D  outline  this 
assumption  in  more  detail. 

5.  The  deviations  from  the  reference  trajectory  are 
small  enough  to  allow  the  use  of  linearized  techniques.  This 
assumption  is  justified  on  the  premise  that  only  small  deviations 
can  be  tolerated  on  this  type  of  mission  without  making  corrective 
techniques  impractical  because  of  the  limited  fuel  capacity  of 
the  vehicle. 

6.  Initial  position  and  velocity  errors  can  be 
predicted  within  a  certain  tolerance.  The  prediction  used  in 
this  study  is  outlined  in  Appendix  D. 

7.  Only  the  outbound  portion  of  a  circumlunar  mission 
between  an  earth  parking  orbit  and  the  moon's  sphere  of  influence 
was  considered.  This  assumption  is  discussed  further  in  Appendix 
A. 

8.  Finally,  the  onboard  computer  is  capable  of  storing 
all  pertinent  constants  and  computing  the  navigation  and  guidance 
problem  in  the  time  allowed. 

Notational  Conventions 

The  notational  conventions  used  in  this  study  are  those  used 
by  Dr.  Battin  (Ref  2:2).  Most  of  the  equations  in  this  report 
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are  three-and  six-dimensional  vector  equations  and  are  denoted 

by  a  capital  letter  such  as  A,  K,  T,  etc.  A  column  matrix  of 

any  dimension  is  represented  by  a  lower  case  underscored  letter 

T 

such  as  a.  A  vector  transpose  is  denoted  ^  and  the  scalar  or 

T 

dot  product  of  two  vectors  a  and  is  written  a  b.  An  overscorce 
denotes  the  average  value,  therefore,  a  is  the  average  value 
of  the  vector  a. 
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I I .  A  Linearized  Navigational  Theory 

A  single  angular  measurement  made  at  a  known  instant  of  time 
is  sufficient  to  fix  the  position  of  the  spacecraft  in  one 
coordinate.  Succeeding  measurements  combined  with  updated  position 
information  from  preceding  measurements  suffice  to  accurately  determine 
the  spacecraft’s  position  and  velocity  (Ref  2:3-4).  The 
practicality  of  this  method  is  based  on  the  fact  that  vehicle 
dynamics  are  governed  by  known  laws  and  the  assumption  that 
deviations  from  a  preplanned  trajectory  are  kept  small  enough  to 
permit  linearization  of  the  problem.  A  discussion  of  the  nature 
of  the  single  angle  fix  and  of  the  linearized  trajectory  equations 
which  constitute  this  method  is  outlined  in  this  section. 

A  Single  Observation  17 ix 

As  stated,  a  single  measurement  serves  to  fix  the  position  of 
the  spacecraft  in  one  coordinate,  if  a  linearized  theory  is 
assumed  (Ref  8:12).  It  is  shown  in  Appendix  B  that  the  deviation 
in  position,  S  of  the  spacecraft  from  the  reference  path  is 
related  to  the  deviation  in  angular  measurement , S An,  by 

^An  =  hnT£rn  (1) 

if  the  observation  is  made  at  a  known  instant  of  time,  tn.  The 
vector  ]ln  is  a  function  of  the  geometry  of  the  celestial  objects 
measured  and  the  type  of  measurement . 
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Because  of  the  inherent  dynamic  coupling  of  position  and 
velocity,  the  result  at  a  later  time,  tn+x,  of  a  measurement 
made  at  tn  does  not  lend  itself  to  a  simple  geometric  interpretation 
However,  to  provide  a  geometric  description  it  is  convenient  to 
define  the  state  of  the  vehicle  dynamics  at  by  the  six  dimensional 
deviation  vector 


S 


Si 

U 


(2) 


where  $  is  the  deviation  in  velocity.  The  relationship 
between  $  an<*  B  is  simply 


^  in+1  “  Tn+l,n^ 


*n 


(3) 


where  ?n+l,  n  is  referred  to  as  the  "state  transition  matrix"  and 
is  defined  later  by  Eq  (19). 

Equation  (I)  may  be  written  in  terms  of  B  as 

SAn  =  hnTKTXxn  (4) 


where  the  rectangular  K  matrix  is  defined  by 


K  = 


6X3 


/ 


By  combining  Eqs.  (3)  and  (4),  we  have 

B  An  *  i?n  K  7n+l,n  B  2n+l 


(5) 


(6) 


Six  observations  made  at  different  times  would  provide  a  set 
of  six  equations  of  the  form  of  Eq.  (6).  If  no  two  of  the 
component  directions  were  parallel  and  the  measuring  device  used 
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to  obtain  S  An  was  perfect,  then  the  deviation  vector  could  be 
accurately  computed. 

Matrix  Solution  of  the  Traj  :tory  Equations 

A  certain  collection  of  so-called  "perturbation"  matrices  are 
basic  to  the  solution  of  the  linearised  trajectory  problem.  The 
following  shows  how  these  matrices  may  be  obtained  as  solutieaa  of 
the  linearized  differential  equations. 

Let  £g(t)  and  v^t)  denote  the  position  and  velocity  of  the 
spacecraft  in  an  earth  centered  inertial  coordinate  system,  and 
let  g(rg,t)  denote  the  gravitational  acceleration.  Then 

dI=  .  Is  .  dI»  (?) 

ar-  5r  * 

are  the  basic  free  fall  equations  of  motion. 

Let  r^/t)  and  v^t)  represent  the  position  and  velocity 
associated  with  the  reference  trajectory  at  time  t  and  define 

r(t)  =  £a(t)  -r^t)  ,  v(t)  =  Tg(t)  -  v^t)  (8) 

Then,  the  deviations  S  r  and  &  v  may  be  approximately  related  by 
the  linearized  differential  equations 

d(  $  £)  =  Si  ,  d(  S  ▼)  _  G(£o,  t)£r  (9) 

dt  dt 

where  G(rc,t)  is  a  3  dimensional  matrix  whose  elements  are  the 
partial  derivatives  of  the  components  of  g'r0}t)  with  respect  to 
the  components  of  r0  ( Ref  15:62). 
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A  fundamental  set  of  solutions  of  Eqs  (9)  nas  been  developed 

(Ref  3:210),  Let  tt  and  t  be  the  time  of  launch  and  arrival, 

L  A 

respectively.  Then  define  the  3  dimensional  matrices  R(t),  R(t)*, 
V(t)  and  V(t)*  as  the  solutions  of  the  matrix  differential 
equations 


dR(  t )  =V(t) 

**<*)«  =  V(t)* 

dt 

dt 

dV(  t )  =  GR(  t ) 

dV(  t )  *  =  GR(  t } 

dt 

dt 

which  satisfy  the  initial  conditions 

R(tL)  =  0  R^V*=  0 

V(tL)  =  I  7'V*=  1 


(11) 


If  we  now  write 


din 

=  R  c 
n  -~ 

+ 

Rn*  £* 

(12) 

n 

=  V  c 
n  — 

+ 

Rn*  £* 

(13) 

where  $ r  =  $ 
— n 

*n 

=  R(tn) ,  etc. , 

and  £  and  £*  are  arbitrary 

constant  vectors,  it  follows  that  these  equations  satisfy  the 
linearized  differential  equations  (9).  They  also  contain  precisely 
the  required  number  of  unspecified  constants  to  meet  any  valid  set 
of  initial  or  boundary  conditions . 

The  elements  of  the  Rn  and  Vn  matrices  represent  deviations 
in  position  and  velocity  along  the  reference  trajectory  resulting 
from  velocity  deviations  at  launch.  A  similar  meaning  is  given 
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to  the  Rn*  and  Vn*  matrices  as  a  resalt  of  deviations  in  velocity 
at  the  arrival  point.  A  complete  description  of  these  perturbation 
matrices  can  be  found  in  the  literature  (Ref  58-66). 

There  are  a  number  of  techniques  for  computing  these 
perturbation  matrices.  Since  this  study  deals  with  an  elliptical 
trajectory,  a  simple  and  straight  forward  matrix  technique 
formulated  by  Dr.  Battin  (Ref.  1:761-773)  was  used.  This 
technique,  which  is  modified  for  the  special  case  of  a  spacecraft 
trajectory  coplanar  to  the  orbit  of  the  moon,  is  outlined  in 
Appendix  C. 

The  State  Transition  Matrix 

The  state  transition  matrix  was  defined  by  Eq  (3)  and  can  be 
computed  from  a  knowledge  of  the  Eqs  (12)  and  (13).  First,  c 
and  c*  must  be  obtained  as  solutions  of  Eqs  (12)  and  (13).  After 
some  simplification 

£  =  -Ln'Vn*  i  r„  4^,)  (14) 

c*  =  -Ln*-l(cn  <f  rn  -  S  »n)  (15) 

where 

Cn  =  Vn  Rn_1 

cn*  =  Vn*  Rn*'1  (16) 

Ln  =  Vn  -  Cn*  Rn 
Ln*  =  Vn*  -  C„  Rn* 

Thus,  with  c  and  c*  determined,  the  position  and  velocity  deviations 
at  any  time  tn  are  given  by  Eqs  (12)  and  (13). 
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In  termst  of  the  six  dimensional  deviation  vector  defined 
by  Eq  (2), 


<f 


X  = 

— n 


-  - 

■ 

Rn  «n* 

c 

v»  Vn*_ 

c* 

(17) 


A  similar  relationship  can  be  written  for  S  -JSn+l*  Considering 
Eq  (3)  we  have 


^  2n+l 


Using  Eqs  (16),  we  have 


Rn+1 

* 

Rn+1 

~Rn 

i - 

* 

c 

cc 

Vi 

v8+i 

vn 

V  * 
vn 

1  -1 


S 


(18) 


n+l,n 


R 


n+1  Rn+1* 


*_  1  -1 

<cn  Ln> 


V„+i  Vn+1* 
Further  simplification  gives 


n+l,n 


0 

-1  -1 
(C„  Ln*) 


-I  C  * 


-1 


-I  C, 


-1 


n 


(19) 


where 


Rn+1  Rn+1# 

"(Rn  -  (An*Vn  -  Rn)'^,,* 

vn+l  Vn+1* 

(Rn*  *  AnO'1  (AnVn*  -  Rn*)  _1An 

(Ref  3:305) 

-I 

-1 

=  Cn  =  RnVn 

(20) 


A  «  s  C  *"1  -  R*\ n*"1 
n  n  n 


(21) 


The  form  of  Eq  (2">)  was  motivated  by  the  values  of  R(tL)  and 
R*(tA)»  From  Eqs  (16)  it  follows  that  these  two  singular  matrices 
would  make  the  computation  of  C(t^)  and  C(tA)#  impossible.  However, 
this  difficulty  is  avoided  by  the  use  of  Eq  (20) . 
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An  interesting  feature  of  the  state  transition  matrix, 

which  allows  the  checking  of  its  computation,  is  that  it  is 

an  example  of  a  symplectic  matrix  (Ref  3:306)  .  An  even-dimensional 

matrix,  A,  is  said  to  be  symplectic  if 

T 

A  J  A  =  J  (22) 


where 


J 


I 

0 


(23) 


In  addition,  the  absolute  value  of  the  determinate  of  a  symplectic 
matrix  is  equal  to  unity.  These  features  were  used  to  check  the 
computation  of  the  state  transition  matrices  computed  in  this 
study . 
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III.  An  Estimate  of  Position  and  Velocity  Deviation 

The  previous  section  described  a  simple  position  and  v  .ocity 
determining  technique  which  made  no  provision  for  the  inherent 
errors  common  to  all  navigational  measuring  instruments.  Usually 
it  is  necessary  to  discuss  navigational  accuracy  on  a  probabilistic 
rather  than  absolute  basis  because  of  the  inexact  knowledge  of 
the  error  coefficients  of  the  various  measuring  devices  in  the 
system.  Therefore,  a  statistical  evaluation  of  the  accuracy  of 
the  navigation  process  must  be  made.  In  this  study  only  the 
optical  measuring  device  was  considered  as  a  source  of  the 
uncertainty.  Accurate  estimates  of  position  and  velocity  along 
a  space  trajectory  can  be  obtained  by  applying  the  method  of  least 
squares  to  the  system's  errors.  A  linear  least  squares  estimation 
procedure,  where  all  statistical  calculations  are  based  on  first 
and  second  order  averages,  is  outlined  in  this  section  (Ref  2:4-6). 

St  > ' •  £  ;ical  Parameters  and  Definitions 

Berore  discussion  of  the  statistical  parameters,  it  is 
necessary  to  define  some  terms 

$ An  =  the  true  deviation  in  the  angle  An  from  its  reference 
value  at  tn 

£  A„  =  the  observed  or  measured  deviation 
n 

r*  a 

a  A  =  the  estimated  deviation 
n 

<f  A^  ~  Sk  extrapolated  from  a  previous  estimate  by  the 
1  state  transition  matrix  (T) 

an  =  the  error  in  the  measurement  of  S An 
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=  the  error  in  position  estimate 
d^  =  the  error  in  velocity  estimate 
Similar  definitions  can  be  applied  to  Eq  (2). 

From  these  definitions  we  have, 

^An  =  ^An  +  an  <24) 

In  the  subsequent  analysis,  an  will  be  assumed  to  be  a  random 
(a^  =  o)  and  independent  ( an-latT  =  °)  variable.  The  application 
of  the  Central  Limit  Theorem  indicates  that  the  probability 
distribution  function  for  the  system  error  will  be  normally  distributed 
(Ref  9:315).  Furthermore,  the  standard  deviation  ( (f  n)  for  the 
system  error  will  be  the  root- sum- square  combination  of  the 
standard  deviation  for  the  individual  error  coefficients.  Therefore, 
in  this  study  the  variance  of  the  optical  measuring  device  can 
be  computed  by 

2  _ _ 

cTn  =  an2  (25) 


In  the  next  sub-section,  an  estimation  procedure  is  discussed 
for  determining  an  optimum  linear  estimate  of  ^  2En»  denoted 
by  J  xn.  An  integral  part  of  this  estimation  technique  is  the 
correlation  matrix  of  the  measurement  errors  in  the  estimate. 

If  we  write 


^  *n  “  $ 


— n 


£n 


(26) 


then 


(27) 
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The  correlation  matrix  of  the  measurement  errors  is  thus  defined  by 


* 

~ 

_ 

_  T 

E  =  e  e  - 

n  — rr~n 

f 

- 

E  1 
n 

t 

T 

d  d 
— n— n 

3 

E 

n 

4 

A 

n 

1  4 

where  the  traces  of  ER  and  En  are  the  wean-squared  errors  in 
the  estimate  of  position  and  velocity,  respectively  or  as  it  is 
used  in  this  report,  the  uncertainties  in  position  and  velocity. 
Again  from  the  definitions  we  have 

SK  *  Vn-lSVi  <29> 


In  like  manner,  the  extrapolated  correlation  matrix  can  be 
shown  to  be 


V  *  Tn,n-1  Bn-1  Tn  ,0-1  <29A> 

Improving  the  Estimate 

The  development  of  an  optimum  linear  estimate  as  a  recursive 
operation  in  which  the  current  best  estimate  is  combined  with 
newly  acquired  information  to  produce  a  better  estimate  was 
presented  by  R.E.  Kalman  (Ref  12)  and  completely  formulated  by 
R.  H.  Battin  (Ref  2:16-18).  This  development  assumes  that  an 
initial  En  ^  (Appendix  D)  is  known  and  that  a  single  navigational 
measurement  of  the  type  described  in  Appendix  B  is  made  at  time 
tQ.  A  linear  estimate  for  the  deviation  vector  S  xn  expressible 
as 


£*n  =  +  ^  An  <30) 

where  is  a  weighting  factor  which  will  be  chosen  so  as  to 

minimize  the  mean-square  error  in  the  estimate.  It  can  be 
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shown  that  the  correlation  matrix,  defined  by  Eq  (28)  muv  be 
expressed  as  a  function  of  w^  as 

EnCSn)  =  Cl  -  Mn^TKT )  En(I  -  Kh^T)  ♦  Un2nTan2  (31) 


In  order  to  determine  the  optimum  w^  the  technique  of 
variation  calculus  can  be  shown  to  yield 

S  °n2Ciin)  =  2tr  -S  W*T*n'  U  -  %£aT)  + 

1^  en2(iin)  is  to  vanish  for  all  variations  c)  *n*  then  it  must 


(32) 


follow  that 


cnSn  =  *n'  %  (33) 

where  the  scalar  coefficient  cn  is  computed  from 

°n  =  0%'  %d  -  «n2  04) 

Having  obtained  the  optimum  w^,  Ejj  may  be  written  as 

Bn  *  -V^Sn'  «n)  C*n'  Bfen)1  05) 

Equations  (30)  and  (35)  then  serve  as  recursive  relations  to  be 
vised  in  obtaining  improved  estimates  of  position  and  velocity 
deviations  at  each  of  the  measurement  times  selected  along  the 
spacecraft’s  trajectory. 
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IV.  Results 

The  effecf  position  sad  velocity  uncertainties  by  use  of 
various  navigation  fix  sehssses,  altering  the  inclination  eagle 
between  the  reference  trajectory  plane  and  the  lunar  plane, 
considering  different  values  for  the  variance  in  the  errors  of  the 
optical  sssesuring  device,  and  using  different  time  intervals  between 
each  navigation  fir  is  outlined  is  this  chapter*  -Finally,  a  comparison 
between  ?iac  1,  vfai oh  uses  ©sly  esrih-asoon  angular  seasureawnts  for 
input  into  the  recursive  navigation  equations  and  ?ix  2,  which  considers 
and  selects  the  best  of  sarth-atar,  soon»stars  sad  earth-moon  angular 
s^mirasaants  is  presented  in  the  last  section  of  this  chapter. 

Since  only  the  outgoing  elliptical  portion  of  the  lunar  trajectory 
is  considered  in  this  study,  the  final  position  and  vaioeity 
uncertainties  as  listed  actually  represent  the  initial  uncertainties 
for  the  selenocentric  hyperbolic  portion  of  the  lunar  trajectory.  A 
value  of  one  statute  sells  (Bsf  5*13)  was  selected  as  the  desired 
mTinaag  uncertainty  in  final  position  and  0*025  *ph  was  chosen  as  the 
desired  m&xiagsm  uncertainty  in  final  velocity.  This  desired  final 
Etaxisusa  velocity  corresponds  to  approximately  0*001*5  of  the  velocity 
of  the  vetilale  at  that  point  (JW  6  ?3I ;  fW  16  s  755 )  *  It  was  also 
desired  to  maintain  the  mxiwst  position  and  velocity  uncertainties 
during  the  flight  within  an  order  of  magnitude  of  those  found  by 
Sat tin  (8a?  2 f 35-36 ) -  To  help  indicate  the  general  shapes  of  the 
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curves,  the  tables  in  this  chapter  also  include  values  at  the  midcourse 
point,  which  for  this  study  is  considered  that  point  corresponding  to 
25  hours  of  flight* 

Effect  of  Navigation  Measurements 

The  effect  of  using  various  navigation  fix  schemes  on  the 
uncertainties  in  position  and  velocity  is  presented  in  Fig*  1  and 
Table  I.  Fix  1  uses  only  earth -moon  angular  measurements  for  input 
into  the  recursive  navigation  equations;  Fix  2  considers  and  selects 
the  best  of  the  earth-star,  moon-star,  and  earth -moon  angular 
measurements;  Fix  3  considers  and  selects  the  best  of  the  apparent 
diameter  measurements  of  the  earth  and  moon,  and  the  earth-moon 
angular  measurements;  and  Fix  4  considers  and  solects  the  bes*  from 
all  of  the  above  measurements.  For  each  fix,  the  beat  'measurement 
is  considered  that  measurement  which  results  in  the  smallest  value 
for  uncertainty  in  position.  The  measurements  are  further  discussed 
in  Appendix  B.  The  inclination  angle  between  the  lunar  plane  and 
the  reference  trajectory  plane  is  900  (Appendix  A)  and  the  time 
interval  between  each  navigation  fix  is  held  constant  at  3°  minutes* 

As  indicated  in  Fig.  1,  the  position  and  velocity  uncertainty  curves 
each  use  different  values  for  the  variance  in  the  errors  of  the 
optical  measuring  device. 

The  position  uncertainty  curve  for  Fix  3  exhibits  unexpected 
results  which  have  not  been  explained,  although  the  pertinent  equations 
and  computer  program  have  been  checkea  for  possible  errors. 
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Table  I 

The  Effect  of  Navigational  Measurements  on  Uncertainty  in  Position 
and  Velocity 

Note:  This  data  is  taken  from  the  same  computer  results  used  to 
plot  Fig.  1. 


Nav . 

Position  Uncertainty 

(mi. ) 

Velocity  Uncertainty 

(mph) 

Fix 

Maximum 

Mid. 

Final 

Maximum 

Mid. 

Final 

1 

13.667  (2.5)a 

2.841 

.951 

9.829  (l.)a 

.1411 

.0189 

2 

2.092  (14. )a 

1.930 

.787 

5.161  ( .5)* 

.0959 

.0152 

3 

9.321  (9.0)a 

3.222 

.972 

b 

b 

b 

4 

2.081  (16. )a 

1.930 

.787 

b 

b 

b 

a  Time  in  hours  after  time  of  launch  when  maximum  occured. 
b  No  data  taken. 


Table  II 

The  Effect  of  the  Inclination  of  the  Lunar  and  Spacecraft  Planes  on 
the  Uncertainty  in  Position  and  Velocity 

Note:  This  data  is  taken  from  the  same  computer  results  used  to 
plot  Fig.  2. 


llnclin . 

'  Angle 
.(degrees) 

A - - - 

Position  Uncertainty  (mi) 

Velocity  Uncertainty  (raph) 

Maximum 

Mid. 

Final 

Maximum 

Mid. 

Final 

1  0 

56.610  ( 14 . 8 ) a 

13.946 

2.331 

15.349  ( .25)a 

.7609 

.0597 

30 

15.264  (1.5)a 

10.677 

4.075 

17.057  ( .25)a 

.5137 

.0539 

60 

39.166  ( 4 . 8) a 

3.450 

1 .306 

15.830  ( . 25 ) 3 

.1809 

.0243 

90 

5.693  ( 4 .3) a 

1 .916 

0 . 665 

16.126  (  .25)a 

.1008 

.0135 

a  See  Table  I . 
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A*  expected  Fix  1  has  the  largest  uncertainty  in  position  (13.667 
ailes)  since  it  uses  only  earth-moon  angular  measurements,  but 
the  curve  for  Fix  3  which  also  considers  apparent  earth-moon 
diameter  measurements  crosses  and  remains  above  the  former  curve 
after  7.5  hours  of  flight.  Because  of  these  results,  the  earth 
and  moon  diameter  measurements  are  not  considered  in  the  remaining 
sections  of  this  chapter. 

All  final  position  uncertainty  values  are  below  one  mile.  The 
velocity  uncertainty  curve  for  Fix  1  remains  above  the  curve  for 
Fix  2  for  all  flight  times  and  the  Fix  1  curve  reaches  a  maximum 
of  9'. 829  mph,  compared  to  a  maximum  of  5.160  aph  for  Fix  2.  Both 
final  velocity  uncertainty  values  fall  below  0.025  mph. 

Effect  of  Inclination 

The  effect  of  the  use  of  coplanar  and  non-coplanar  trajectories 
on  position  and  velocity  uncertainties  using  Fix  1  is  shown  in  Fig. 

2  and  Table  II.  Inclination  angles  of  0°,  30°,  60°,  and  90°  are 
considered.  The  time  spacing  between  each  navigation  fix  is  held 
coastant  at  15  minutes  and  a  constant  variance  in  the  errors  of  the 
measuring  device  is  assumed.  Curves  for  Fix  2  and  a  non-coplanar 
reference  trajectory  are  also  shown  for  comparison  purposes. 

The  coplanar  trajectory  resulted  in  the  highest  maximum  value 
for  the  uncertainty  in  position  (56.610  miles)  but  the  curve  drops 
below  the  curve  corresponding  to  an  inclination  angle  of  30°  after 
30.25  hours  of  flight.  For  an  inclination  angle  of  90°,  the  position 
uncertainty  curve  reaches  a  maximum  value  of  only  5.693  miles  and  it 
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.i  *  60° 


gl  'i  *  90°Y 
1  \  x 


0  X©  20  30  40  ‘ 

Plight  Tine  (hears) 

(b)  If  fact  an  Veleeity  Uncertainty 

Pig  2.  The  Iff  act  of  Inclination  af  the  Lunar  and  fihference 
Trajectory  Planes  an  Uncertainty  in  Position  and 
Velocity.  (Pin  l.Atf*  1$  sin,  (Xn2-  (O.OOOQ5)2) 
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is  the  only  curve  to  drop  below  one  mile.  The  maximum  velocity 
uncertainties  compare  in  order  of  magnitude  for  all  values  of 
inclination  angles  considered  but  the  curve  for  90°  drops  much 
faster  than  the  others.  The  curve  for  all  inclination  angle  of 
30°  crosses  above  the  curve  for  60°  after  14.50  hours  of  flight. 

Only  the  final  velocity  uncertainties  for  inclination  angles  of 
60°  and  90°  are  below  0.025  mph. 

Effect  of  Fix  Spacing 

The  effect  on  position  and  velocity  uncertainties  by  using 
different  time  spacings  between  each  navigation  fix  is  shown  in 
Fig.  3  and  Table  III.  These  values  are  for  Fix  1  using  an  inclination 
angle  of  90°  and  assuming  a  constant  variance  in  the  errors  of  the 
optical  measuring  device.  Curves  and  data  are  shown  for  fix  time 
intervals  of  one  hour,  30  minutes,  15  minutes,  5  minutes,  and  one 
minute . 

From  the  curves  it  can  be  seen  that  the  shorter  the  time  spacing 
between  each  navigation  fix,  the  lower  the  uncertainties  in  both 
position  and  velocity.  For  a  time  spacing  between  each  fix  of  one 
minute,  the  maximum  uncertainty  in  position  is  only  0.430  miles 
as  compared  to  12.34  miles  for  a  one  hour  time  spacing.  Only  the 
curves  for  time  intervals  of  one  and  five  minutes  fall  below  a  final 
value  of  one  mile.  In  general  the  velocity  uncertainty  curves  follow 
the  same  pattern  as  the  position  curves ,  although  their  maximum  values 
do  not.  The  maximum  velocity  uncertainty  of  16.126  mph  occurs  on 
the  curve  for  a  fix  spacing  of  15  minutes  and  the  smallest  maximum  of 
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Table  III 


The  Effect  of  Fix  Spacing  on  the  Uncertainty  in  Position  and  Velocity 

Note:  This  data  is  taken  from  the  same  computer  results  used  to 
plot  Fig.  3. 


A  tf 

(Hours) 

Position  Uncertainty  (mi.) 

Velocity  Uncertainty  (mph) 

Maximum 

Final 

Maximum 

Mid. 

Final 

1 

12.337  (3.0)a 

3.709 

1.295 

9.892  (1.0)a 

.1946 

.0260 

h 

11.392  ( 2 .0) a 

2.686 

.932 

9.829  (1.0)a 

.1411 

.0189 

k 

5.693  (4.3)a 

1.916 

.665 

16.126  ( .25)* 

.1008 

.0135 

1/12 

1.681  (  .17)a 

0.864 

.366 

12.527  (  .17)a 

.0454 

.0075 

|  1/60 

0.430  ( 4 .4) a 

0.203 

.054 

11.363  ( .Q2)a 

.0109 

.0014 

a  See  Table  I 


Table  IV 

The  Effect  of  Sextant  Error  on  the  Uncertainty  in  Position  and  Velocity 

Note:  This  data  is  taken  from  the  same  computer  results  used  to 
plot  Fig.  4. 


Sextant 

Position  Uncertainty  (mi.) 

Velocity  Uncertainty  (mph) 

Errorb 

Maximum 

Mid. 

Maximum 

Mid. 

Final 

CTna2 

9.304  ( 2 .8) a 

2.032 

.678 

16.430  (  .25)a 

.1067 

.0137 

<r  nb2 

6.873  (2.8)a 

1.437 

.480 

16.272  C .25)a 

.0755 

.0097 

6  nc 

5.693  ( 4.3) a 

1.916 

.665 

16.126  (  .25) a 

.1008 

.0135 

.  2 
<5  nd 

4.050  (4.0)a 

1.355 

.470 

16.095  ( .25)a 

.0713 

.0096 

a  See  Table  I 

p  2  2 

b  {Tna  is  defined  by  Eq.  D-2  in  Appendix  D,  O'njj  =%  <f na 
<fnc2  =(0.00005)2  rad.,  (f  nc2 
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9.829  mph  occurs  on  the  curve  for  a  spacing  of  30  minutes.  All 
final  velocity  uncertainties  are  below  0.025  mph  ex'*»pt  for  the 
curve  for  a  one  hour  fix  spacing. 

Effect  of  Sextant  Errors 

The  effect  of  using  different  values  for  the  variance  in  the 
errors  of  the  optical  measuring  device  on  uncertainty  in  position 
and  velocity  is  shown  in  Fig.  4  and  Table  IV  using  Fix  1.  The  time 
interval  between  each  fix  is  held  constant  at  15  minutes  and  the 
reference  trajectory  plane  is  inclined  to  the  lunar  plane  at  90°. 

Four  different  values  for  the  variance  are  considered  (Table  IV, 

Note  b) . 

As  expected,  the  maximum  uncertainty  in  position  of  9.304 

2 

miles  occurs  on  the  curve  related  to  {f  n  ,  as  this  curve  accounts 

for  the  larger  uncertainty  in  defining  the  horizon  when  the  vehicle 

is  near  the  earth  or  moon.  The  maximum  position  uncertainty  for  the 

2 

constant  variance  <f  nc  is  5.693  miles.  After  about  25  hours  of 

2  2 

flight,  the  curves  corresponding  to  (f  na  and  cf  nc  are  almost 

2  2 

identical  as  are  the  curves  for  <S  n^  and  cf  nc  .  All  position 

uncertainties  at  arrival  are  less  than  one  mile  and  the  velocity 

uncertainty  curves  follow  a  pattern  similar  to  the  position  curves. 

2 

The  maximum  velocity  uncertainty  occurs  on  the  curve  for  Cf  na  at 

.  2 

16.430  mph  and  the  minimum  is  on  the  curve  for  0  n^  at  a  value  of 
16.095  mph.  All  final  velocity  uncertainties  are  less  than  0.025  mph. 
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(a)  Effact  an  Paaitlaa  Uncartainty 


(b)  Zffaet  on  Talaoity  Uncartainty 

4*  Xff aat  of  Saoctamt  Irrara  an  tka  Uncartainty  in 
Paaitian  and  Yalacity.  (fix  1,  i*90*,At^*  miu) 
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Comparison  of  Fix  1  and  Fix  2 

Fig.  5  and  Table  V  show  a  coplanar  comparison  between  Fix  1 

and  Fix  2.  Fix  2,  with  a  time  spacing  between  each  measurement 

held  constant  at  15  minutes,  is  considered  the  standard.  Each  fix 
2 

uses  (S  na  as  the  variance  in  the  errors  of  the  optical  measuring 
device.  Fix  1  is  plotted  for  constant  fix  spacings  of  one  minute 
and  30  seconds  respectively.  The  former  curve  has  a  maximum  position 
uncertainty  of  14.634  miles  and  the  latter  10.095  miles  as  compared 
to  7.011  miles  for  Fix  2.  Both  position  uncertainty  curves  for  Fix 
1  drop  below  the  curve  for  Fix  2  before  40  hours  of  flight  and  all 
final  position  uncertainties  are  less  than  one  mile.  Fix  1,  with  a 
time  spacing  between  each  fix  of  1  minute,  has  a  maximum  velocity 
uncertainty  of  12.878  raph  and  the  same  Fix,  with  a  time  spacing 
between  each  fix  of  ^0  seconds,  has  a  maximum  velocity  uncertainty  of 
12.614  mph  as  compared  to  14.953  mph  for  Fix  2.  Both  velocity 
uncertainty  curves  for  Fix  1  exhibit  a  plateau  region  between  times  of 
flight  of  approximately  3  to  15  hours  which  are  probably  the  result 
of  the  non-optimum  nature  of  the  earth-moon  angular  measurements  during 
the  early  part  of  the  flight.  All  final  velocity  uncertainties 
axe  less  than  0.025  mph. 

Fig.  6  and  Table  VI  show  a  non-coplanar  comparison  between 
Fix  1  and  Fix  2.  The  inclination  angle  is  90°  and  the  same  variance 
in  the  errors  of  the  optical  measuring  device  is  assumed  as  used 
above.  The  time  spacing  for  Fix  2  is  again  15  minutes  and  Fix  1  is 
plotted  for  time  spacings  of  five  and  one  minutes.  The  maximum 
position  uncertainty  for  Fix  1,  with  a  time  spacing  between  measurements 
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(a)  Iffect  on  Position  Uncertainty 


(b)  Effect  on  Velocity  Uncertainty 
fig,  5,  Co  planar  Coapariaon  of  Fix  1  and  Fix  2*  (  i*0°# 

O' 
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Table  V 

Coplanar  Comparison  of  Fix  1  and  Fix  2  (i~0°) 


Note:  This  data  .la  taken  from  t.  ,  same  computer  results  uses  to 
plot  Fig.  3* 


N a v. 

Posi.ion  Uncertainty  (ad)  j 

; 

,  Velocity  Uncertainty  (mph)  j 

Fix 

Maximum 

Mid 

Final 

Final  | 

lb 

14.634  (13.8)a 

3.735 

0.620 

12.878  (0.l8)a 

— 

0.2038 

0.0139* 

!d 

10.093  (11.9)a 

2.640 

0*438 

12.614  (0.13)a 

0.1440 

0.0112 

2C 

7.011  Q.?3)a 

2.333 

0.738 

14.953  (o.23)a 

0.1260 

0.0200 

a  See  Table  I  c  At^=  13  min 

b  A  I?-  1  min  d  A  j)0  sr  c 


Table  VI 

Non-coplarar  Comparison  of  Fix  1  and  Fix  2  (i  =  90°) 


Note:  This  df; 'a  is  taknn  from  the  same  computer  roe: 'Its  used  to 
plot  Fig.  6. 


r 

Nav, 

'1 

Position  Uncertainty  (mi) 

—  *  ■  -■  L„  ...  .... - — - — 

Velocity  Uncertainty  (>apb) 

Fix 

! 

Maximum  i 

Mid 

Final 

Maximum 

Mid 

Final 

lb 

- j. 

2.409  (7t08)a  | 

I 

1.148 

0.392 

12.633  (o.i?)a 

0.0604 

Q.0080 

1.302  (0.32)a  ! 

0.216 

0.079 

11.364  (0e02)a 

0.0116 

0,0018 

L_2° _ ! 

1.507  (12.7)k  j 

1.382 

0.339 

!  6.923  (0.23)a 
_ 

0 ,0693 

L 

0.0108 

J 

a  See  Table  I  e  At^s  13  min 

b  A  t^=  5  min  ci  At  -  1  min 

1  f 
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of  five  minutos,  is  2.409  miles  and  the  maximum  for  the  same  Fix, 
with  a  time  spacing  of  one  minute,  is  1.302  miles  as  compared  to 
1.507  miles  for  Fix  2.  All  final  position  uncertainties  are  less 
than  one  mile.  The  maximum  velocity  uncertainties  for  Fix  1  are 
12,655  raph  and  11.3638  mph  as  compared  to  6.9226  mph  for  Fix  2. 
The  velocity  uncertainty  curve  for  Fix  1,  with  a  time  spacing  of 
five  minutes,  shows  a  plateau  region  early  in  the  flight  as  was 
found  with  the  coplanar  curves.  All  final  velocity  uncertainties 
are  less  than  0.025  mph. 
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V.  Conclusions 

The  requirement  for  estimating  the  position  and  velocity 
accurately  along  an  earth-moon  freefall  trajectory  using  only  earth 
moon  angular  measurements  and  an  on-board  computer  capable  of 
computing  the  recursive  navigation  equations  has  been  met.  As 
shown  in  the  preceding  section,  updating  the  navigation  information 
every  minute  or  less  for  the  coplanar  reference  trajectory  and 
every  five  minutes  or  less  for  the  non-coplanar  reference 
trajectory  will  result  in  acceptable  uncertainties  in  the  estimate 
of  position  and  velocity. 

With  the  use  of  non-coplanar  reference  trajectories,  a  tra'  - 
off  should  be  obtainable  between  the  capability  of  the  on-board 
computer  and  the  desired  uncertainties  in  position  and  velocity 
measurements.  This  study  has  shown  that  the  use  of  only  earth- 
moon  angular  measurements  with  an  on-board  recursive  navigation 
system  is  a  feasible  navigation  technique  for  a  lunar  mission. 
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VI .  Recommendations 

Various  extensions  to  this  study  are  suggested  by  Dr.  Battin's 

work  as  outlined  in  iterances  2  and  3.  However,  the  next  logical 

development  of  this  study  would  be  the  estimation  of  velocity 

correction  along  the  flight  path  needed  to  reach  the  target  point 

(Ref  2:21-36).  In  this  extended  study,  full  use  could  be  made  of 

the  perturbation  matrices,  state  transition  matrices,  and  position 

and  velocity  estimates  generated  for  this  study.  The  digital 

programs  used  in  this  study  can  be  obtained  from  Mr.  Poole  of  the 

A.S.D.  Digital  Computation  Division  (see  Preface)  under  problem 

number  64024  and  task  number  040A-77137.  It  is  suggested  that  any 

new  studies  be  based  on  the  data  forAtf=  5  min.  or  less,  inclination  = 

2  2 

60°  or  90°,  and  6  n  =  (0.00005)  . 


37 


nnr/ittt  /64.-Q 


Bibliography 


1.  Battin,  R.H.  A  Recoverable  Interplanetary  Space  Probe , 

Volume  IV.  Report  R-235.  Cambridge,  Mass:  M.I.T. 

Instrumentation  Laboratory,  July  1959 

2.  Battin,  R.H.  A  Statistical  Optimizing  Navigation  Procedure 
for  Space  Flight .  Report  R-341.  Cambridge  Mass:  M.I.T. 
Instrumentation  Laboratory,  September  1961  (Revised  May  1962). 

3.  Battin,  R.H.  Astronautical  Guidance.  New  York:  McGraw-Hill 
Book  Co.,  1964. 

4.  Buchheim,  R.W.  Motion  of  a  Small  Body  in  Earth-Moon  Space . 

Rand  Research  Memo.,  RM-1726.  Santa  Monica,  California: 

The  Rand  Corporation,  June  1956. 

5.  Danby,  J.M.A.  "Matrix  Methods  in  the  Calculation  and  Analysis 
of  Orbits."  A. I. A. A.  Journal,  2:  13-16  (January  1964). 

6.  Duncan.  Donald  B.  "Inertial  Guidance  and  Space  Navigation." 
Journal  of  the  Institute  of  Navigation,  6z  30-33,  (Spring  1958) . 

7.  Frimtzis,  R.  Lunar  Vehicle  Guidance  Study .  Technical 
Documentary  Report  ASD-TDR-62-2D7  San  Diego,  California 
General  Dynamics/Astronautics f  March  1962. 

8.  Hamer  H.A.  and  A.P.  Mayo.  Error  Analysis  of  Several  Methods 
of  Determining  Vehicle  Position  in  Earth-Moon  Space  from 
Simultaneous  Onboard  Optical  Measurements .  NASA  Report 

TN  D-1805 .  Hampton,  Va:  Langley  Research  Center,  June  1963. 

9.  Harter,  G.A.  Error  Analysis  and  Performance  Optimization  of 
Rocket  Vehicle  Guidance  Systems ,  Chapter  12,  Inertial  Guidance 
Edited  by  G.R.  Pitman.  New  York:  John  Wiley  and  Sons  lac.,  1962. 

10.  Huss,  C.R.  H.A.  Hamer s  and  J.P.  Mayer.  Parameter  Study  ox 
Insertion  Conditions  for  Lunar  Missions  Including  Various 
Trajectory  Considerations .  NASA  Report  TR  R-122.  Langley 
AFB,.Va:  Langley  Research  Center,  1961. 

11.  Jones,  R.L.  and  A.P.  Mayo.  A  Study  of  Some  Transition 
Matrix  Assumptions  in  Circumlunar  Navigation  Theory.  NASA 
Report  TN  D-1812.  Hampton,  Va:  Langley  Research  Center, 

October  1963. 

12.  Kalmar.,  R.E.  "A  New  Approach  to  Linear  Filtering  and 
Prediction  Problems."  Transactions  of  the  A.S.M.E. ,  Series 
D.  82:  35-45,  (March  1900). 


GQC/EE/64-9 


13.  Michael,  W.H,  Jr.  and  R.H.  Tolson.  Effect  of  Eccentric ' ty  of 
the  Lunar  Orbit ,  Obi ateness  of  the  Earth ,  and  Solar  Gravita¬ 
tional  Field  on  Lune  r  Trajectories .  NASA  Report  TN  D-2?7  . 

Langley  Field,  Va:  Langley  Research  Center,  June  1960. 

14.  Schwaniger,  A  J.  Trajectories  in  the  Earth-Moon  Space  with 
Symmetrical  Free  Return  Properties*  NASA  Report  TN  D-1833. 
Huntsville ,  Ala:  George  C.  Marshall  Space  Flight  Center, 

June  1963. 

15.  Scott,  D.R.,C.J.  Januska,  and  R.E.  Willes.  Optimum  Statistical 
Operations  with  Celestial  Fix  Data  for  Interplanetar y  Navigation . 
Unpublished  M.S.  Thesis  (T-298) .  Massachusetts  Institute  of 
Technology,  June  1962. 

16.  Swanlund,  George,  "Analysis  Techniques  to  Determine  Guidance 
Computation  Requirements  for  Space  Vehicles."  American 
Rocket  Society  Journal,  32:  755-761  (May  1962). 


39 


GGC/EE/ 64-9 


Appendix  A 

Lunar  Reference  Trajectory 

The  utilization  of  trajectory  shapes  which  bring  an  exploring 
vehicle  from  the  earth  to  a  point  arbritrarily  near  the  moon  and, 
which  if  allowed  to  continue,  return  the  lunar  vehicle  to  the 
earth  with  such  conditions  that  reentry  and  recovery  are  feasible, 
might  be  attractive  for  both  instrumented  and  manned  flights. 
Photographing  the  surface  of  the  moon,  determination  of  astronomical 
constants  by  observation  of  trajectory  perturbation,  and  checkout 
of  hardware  are  examples  of  unmanned  lunar  missions.  The  first 
manned  exploration  flights  may  be  "free-return"  circumlunar  flights 
with  no  plan  to  land  on  the  moon.  If  the  manned  mission  is  to 
land  on  the  moon,  a  ''free-return*’  trajectory  nay  be  used  so  that 
if  difficulties  arise  which  would  make  the  lunar  landing  undesirable 
the  astronauts  will  return  safely  to  the  earth  (Ref  14;l-2) . 

Because  this  is  a  feasibility  study  and  to  simplify  the 
reference  trajectory  calculations,  several  simplifying  assumptions 
relating  to  the  mathematical  model  of  the  earth-moon  system  have 
been  made. 

(1)  The  earth  and  moon  are  homogenous  spheres  and  are 
considered  point  masses . 

(2)  The  earth  is  located  at  the  center  of  its  sphere  of 
influence  and  when  the  lunar  vehicle  is  traveling  inside  this  sphere, 
its  ”f ree-return"  trajectory  is  effected  only  by  the  gravitational 
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field  of  the  earth.  Thus  during  this  portion  of  flight,  the 
path  of  the  spacecraft  describes  a  conic  with  the  earth  at  one 
focus  which  will  be  elliptic,  parabolic,  or  hyperbolic  depending 
on  the  insertion  parameters.  An  analogous  situation  exists  after 
the  spaceship  enters  the  sphere  of  influence  of  the  moon  with  the 
exception  that  the  path  of  the  vehicle  will  always  be  hyperbolic. 

The  approximate  radius  of  the  sphere  of  influence  of  the  moon  has 
been  calculated  to  be  35,781  miles  (Ref  10:29). 

(3)  The  reference  trajectory  is  high  enough  above  the 
earth  that  the  effects  of  the  atmospheric  drag  may  be  neglected. 

(4)  The  moon  revolves  about  the  center  of  the  earth  in 
the  lunar  plane  with  a  circular  orbit  of  238,857  statute  miles 
(one  lunar  unit)  and  at  an  angular  rate  of  0.22997084  rad/ day 
(Ref  10:29) . 

(5)  The  effects  of  the  sun,  planets,  meteroids,  and 
other  perturbing  bodies  are  neglected. 

(6)  The  effect  of  the  pressure  of  solar  radiation  is 

neglected. 

The  effects  of  the  neglected  factors  on  the  lunar  reference 
trajectory  are  discussed  in  the  literature  and  are  not  presented 
in  this  paper  (Ref  4:61-74;  Ref  11:10-11;  Ref  13).  Deviations 
resulting  from  the  introduction  of  the  neglected  factors  can,  for 
the  most  part,  be  compensated  for  by  changes  in  the  injection 
velocity  (Ref  7:  3-5). 

The  reference  trajectory  can  be  calculated  from  a  knowledge 
of  the  insertion  parameters.  To  increase  the  operational  flexibility 
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of  the  mission,  launch  with  a  non-zero  flight  path  angle  (the 
flight  path  angle  is  defined  as  the  angle  between  the  local 
horizontal  and  the  velocity  vector)  from  a  parking  orbit  will 
be  considered.  The  variable  insertion  conditions  are  the  radius  of 
the  parking  orbit,  flight-path  angle,  magnitude  of  the  insertion 
velocity,  and  the  lunar  lead  angle.  For  a  4100  mile  parking  orbit 
radius,  insertion  velocities  below  the  parabolic  escape  velocity  of 
36,063  ft/sec  will  result  in  elliptical  outgoing  trajectories.  The 
insertion  parameters  selected  for  this  study  (Table  A-l)  will  result 
in  a  coplanar  circumlunar  mission  that  wi  11  bring  the  vehicle  to 
within  1000  miles  of  the  surface  of  the  moon  and  return  the  spaceship 
to  within  40  miles  of  the  surface  of  the  earth  (Ref  8:23) .  Only  the 
outgoing  elliptical  portion  of  the  reference  trajectory  from  launch 
to  the  moon's  sphere  of  influence  is  used  in  this  study  (Fig  A-l) 
as  this  portion  represents  a  large  percent  of  the  total  distance 
traveled  on  a  lunar  mission .  The  equations  for  calculation  of  the 
elliptical  parameters  are  listed  in  Reference  10:  Appendix  B. 

Additional  mission  flexibility  can  be  obtained  if  non-coplanar 
reference  trajectories  are  also  considered.  The  lunar  plane  (x,  y) 
is  inclined  to  the  reference  trajectory  plane  (a,b)  at  various 
inclination  angles  to  obtain  the  effect  of  using  non-coplanar 
reference  trajectories.  This  method  of  inclining  the  two  planes 
eliminates  the  necessity  to  recompute  the  state  transition  matrices. 
Transfovraation  from  the  lunar  coordinates  (x,y.O)  to  the  inertial 
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Trajectory  Insertion  Parameters 


Ineertion 

Parameter 

Symbol 

Talue 

Remarks 

Radius  of 

The  radius  of  the 
intermediate  or 
launch  orbit  as 
measured  from  the 
center  of  the  earth 

Parking  Orbit 

4100  miles 

Plight-path 

angle 

40  degrees 

Angle  between  the 
local  horizon  and 
the  velocity  rector 

Telocity 

35.900 
ft /sec 

Value  chosen  is 

163  ft /sec  less 
than  parabolic 
escape  relocity 

Lunar  lead 

57  degrees 

Angle  between  the 

angle 

X 

radius  rectors  from 
the  center  of  the 
earth  to  the  moon 
and  to  the  launch 
point  at  time  of 
launch. 
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coordinates  (a,b,c),  is  given  by  (The  Inertial  coordinate  system  is 
defined  in  Appendix  C.) 


*  — 

-  — 

a 

X 

b 

SHi 

y 

e 

m  — 

0 

where 


fi 


0 


co  si 


0  siai 


0 

-siai 

cosi 


and  i  is  the  inclination  angle  between  the  lunar  plane  and  the 
reference  trajectory  plane* 
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Appendix  B 

Navigational  Measurements 


The  best  choice  of  celestial  measurements  made  at  any  instant 
of  time  along  a  space  trajectory  depends  on  the  position  of  the 
spacecraft  within  the  geometry  of  the  Space  System*  In  this  study 
a  variety  of  measurements  made  on  different  sets  of  celestial 
objects  in  an  Earth-Moon  System  is  possible.  The  problem  is  to 
define  the  possible  measurements  at  each  instant  of  time  and 
determine  which  measurement  yields  the  best  estimate  in  position 
as  defined  by  Eq  (34)  in  Section  III. 

It  is  assumed  in  this  analysis  that  the  spacecraft  clock  is 
perfect.  This  assumption  is  consistent  with  the  accuracy  of  present 
day  spacecraft  clocks  on  comparatively  short  space  trajectories,  such 
as  the  50  hour  lunar  trajectory  considered  in  this  study.  In 
addition,  the  time  interval  required  to  make  a  single  celestial 
measurement  can  be  neglected,  since  the  measuring  device  can  be 
programmed  to  search  out  and  hold  the  selected  celestial  objects 
prior  to  the  time  of  measurement.  The  relationship  between  position 
error  and  measurement  error  is  essentially  linear  over  a  relatively 
wide  range.  Therefore,  a  linear  perturbation  theory  can  be  used  for 
detailed  studies  of  any  of  the  methods  discussed  in  this  appendix 
(Ref  3:12) . 

It  was  stated  ir?  Section  II  that  each  measurement  establishes 
-  component  of  spacecraft  position  along  some  direction  in  space. 

If  Q  is  the  quantity  to  be  measured  and  $  Q  is  the  difference  between 
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the  true  and  the  reference  values,  then  it  car  be  shown  that  the 
deviation  in  spacecraft  position  i  r  is 

S  Q  =  hT  l  r  (B-l) 

regardless  of  the  type  of  measurement  (Ref  3:223-225).  Thus,  the 
h  vector  alone  will  characterize  the  type  of  measurement. 

In  the  interest  of  computing  simplicity,  the  celestial  objects 
used  in  this  study  were  the  earth,  moon,  and  six  (6)  stars  arbitrarily 
located  at  an  infinite  distance  along  the  coordinate  axes  of  the  space 
coordinate  system. 

Earth-Moon  Measurement 

It,  can  be  shown  (Ref  3:223-225)  that  for  the  earth-moon  measure¬ 
ment 

ill  -(£1-B2)i!2  B2  -  (Hi *H2 >iH 

=  - - - —  +  -»■-  -  ■  —  1  ( B-2 ) 

r  sin  A  2  sin  A 

where  n-^  and  ng  are  the  unit  vectors  from  a  point  on  the  reference 
trajectory  to  the  earth  and  moon  centers  respectively,  r  and  2  are 
the  distances  from  the  reference  point  to  the  earth  anu  moon  centers 
respectively,  and  A  is  the  angle  from  the  earth  line  to  the  moon 
line . 

Planet-Star  Measurement 

This  is  simply  a  special  case  of  the  earth-moon  measurement. 

Eq  B-2  can  be  modified  to  express  the  geometry  vector  of  the  earth- 
star  measurement  as  follows: 
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&S  -fei  *  Ofl) 
r  sin  Ag 


k  - 


(B-3) 


where  rig  is  the  unit  rector  from  a  point  on  the  reference  trajectory 
to  the  star  and  Ag  is  the  angle  from  tls  earth  line  to  the  star 
line. 

Likewise,  a  similar  relationship  can  be  formulated  for  a 
moon-star  measurements 


k8  ~(l2  •  fi2 

£  =  - — - ^ - - 

2  ©in  Aj^ 

where,  is  the  angle  from  the  moon  line  to  the  etar  line. 


(B-4) 


Apparent  Diameter  of  a  Plan&t  Meesureatant 

The  geometry  rector  for  the  apparent  diameter  of  a  planet 
(Ftef  3 5 223-225)  can  be  computed  by 


k  = 


2 

x(4  -  dp2) l/z 


(B-5) 


where  dp  is  the  actual  diameter  of  the  pic  at  (dg  *  4-i»78127X10^  ft. 
and  «  11.40575^0^  ft),  2  is  the  distance  from  the  point  on  the 
reference  trajectory  to  the  planet  center  (r  for  the  earth  and  Z 
for  the  moon) ,  and  is  the  unit  vector  from  reference  point  to 
tha  planet  center  (n^  for  the  earth  and  for  the  moon). 
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Appendix  C 

The  Fundamental  Perturbation  Matrices 
The  procedure  used  to  generate  the  perturbation  matrices 
discussed  in  Section  II  is  developed  in  outline  form  in  Reference 
1:  pgs.  761-773.  Since  the  perturbation  matrices  in  this  study 
were  developed  for  the  elliptical  spacecraft  reference  trajectory 
and  lunar  orbit  in  the  same  plane,  simplifications  to  the  equation 
given  in  the  cited  reference  were  made.  These  simplifications  are 
based  on  the  fact  that  in  this  coplanar  case, 

fl  =  -w  (C-l) 

and  i  =0  (C-2) 

As  shown  in  Fig  C-l,  jfX  is  the  true  anomaly  of  the  spacecraft 
reference  trajectory  at  the  time  of  launch  and  i  is  the  inclination 
of  the  spacecraft  reference  trajectory  and  lunar  orbit  planes 
measured  at  the  ascending  node  (AN)  .  The  parameter  w  is  defined 
for  the  coordinate  systems  used  in  the  cited  reference  and  has  no 
meaning  in  the  coordinate  system  of  Fig.  1C.  It  is  only  stated  here 
as  a  key  to  the  simplication  process. 

Computation  of  R( t )  and  V( t )  Matrices 

The  equations  used  to  compute  the  K{t)  and  V(t)  matrices 
of  this  study  are  lised  here  for  reference.  The  parameters  used 
in  these  equations  are  defined  in  Table  C-I.  The  R(t)  and  V(t) 


matrices  were  written  in  the  form 


CGC/EE/64-9 
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Table  C-I 

Reference  Trajectory  Parameters 


Parameter ( s) 


Definition^ 


U 

t 

tL 

r(L) 

r 

v 

a 

P 

e 

E 

n 

m 


time  of  flight  =  t  -  tL 
time  measured  from  perigee 

time  of  launch  measured  from  perigee  (0.0065147892  days) 
position  of  spacecraft  at  time  t ^  (see  Fig  C-l) 
position  of  spacecraft  at  time  t  (see  Fig  0-1) 
velocity  of  spacecraft  at  time  t 

( 2 \ 

semi-major  axis  (0.95190686  lunar  units'  ') 
semi-latus  rectum  (0.019964129  lunar  units'  1 ) 
eccentricity  (0.98945804) 

eccentric  anomaly  at  time  tL  (0.12312378  rad.) 
true  anomaly  at  time  (1.4052372  rad) 

3  2 

gravitational  parameter  (0.05386804  l.u.  /day  ) 


Note  ^ 
(2) 


The  values  listed  are  computed  from  equations 
obtained  from  Ref  i. 

1  lunar  unit  =  238,857  miles 
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where 


ra(L)  *b(L) 


MppT= 


R  (t)  = 


r(L) 

r(L) 

-rb(L) 

ra(L) 

r(L) 

r(L) 

0 

0 

R11 

r21 

R12 

R22 

(C-4) 


V  (t)  =  V, 


(C- 5) 


The  elements  of  the  R^t)  and  V^t)  matrices  are  as  follows 


R11  =  *  L,  r  .  3  1  a 


WiTimf  3  [r* '  2  v»*f  J  -  a--p  - 

va  |2jr(L)2  -  P(P-r(L))l  ]  +  ^P)*sra(l-)rb 


y«.o  r(L) 


(C-6) 


{yup)h  J  2ae  r  - 
r(L)1 


3  va  tf 
2 


rb2  +  P  a  -  r( L)2 
T~  P 
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va 


i/j ca}%  sin  E 


2ar(L)' 


ZiEz£Lkl)  a 

l 


L. 


rb(L)  rb  (p  +  r(L)) 
(yU,p)h  e  r(L) 


(C-7) 


R12  = 


'22 


a  rb<L> 

^P)V  r(L) 

vb  _ 

tyx. a)^?  sin  E 

=  (AP)^ 

Ae  r(L) 


2e 


rb- 1  v*  *t 


rarb 


0  .  . 2  P(P-r(L)) 

2r(L)  -  - r— 


(C-8) 


e  r(L) 


r 

u 

/  2ae 

_b  - 1  %  ‘f 

- 

t" 

J 

a  -  ^(l)2“ 
P 

l 

»  mm 

L  J 

- 

vb 
sin  e| 


par( L)  -  P( p-r( L) )  a 
e2 


rXL^ 
p  /  !  j 


/  j 


R33  “  Wp) 


rarb^L)  (p  +  r(L)) 
e  (/c?)h  r(L) 


(C-9) 


ra(L)rb  -  ra  rb(L) 


(C-10) 


a  r.(L) 


11  ^p)*£  r(L) 


i*“ 

P 

r 

ra 

3’ 

-  va 

2 

+  va  ra  p 
~"r2 


-  rb  vb  ♦ 


(C-ll) 
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M  xJl. - -  2r(L)' 

r  i/i a)^5  sin  E 


V?-,  =  /^er(L) 


+  Uv)*8  ra;L)vb 

^er(L) 


L  2ae  3/<  £^tf  -  va  +  a  "  iitlf. 

I  O  S  *._/  ->  r-> 


vavaP  "  rbvb 


2 y  r-  2  ! 


e  r3(/ca)^  sin  E 

1 


2ar(L) 


P  (p  -  r(L) )  ( a  -  r(L)2)  t  -  rb(L)vb  (P  +  r(L))  (C-12) 


(/tp)*  e  r  (L) 


a  rb(L) 
^<P)^  r(L) 


-  rb  t,  -  Vk 

2  73  f  2^ 


... ,  b  -  2r(L)' 

( sxn  E 


^a£bE  +  £aYh.  + 

r**  r 


UP)^a(L)va 

( L) 


(C-13) 


V22  ™  yftP.1 — —  V  2ae 

2  ^er(L)  ] 


va  rbP  +  *b  vb 


rb  tf  -v6  + 

2 2" 


a  -  rjLjl 
P 


er3  (^.a)%  sin  E 


2ar(L)2  -  (p  _  ;,(L})  (a  -  rjLj2  *  rb(L)va(p  +  r(L)) 
C  P  -  e  ( /<  p) ^  r ( L) 


r 

V33  =  jX’a{L)Vb  **  Va  rb(L) 


(C-14 

(C-15) 
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Computation  of  B»(t)  and  V»(t)  Matrices 

The  R*(t)  and  Y*(t)  matrices  are  computed  using  the  same  equations 
derived  for  R(t)  and  V(t)  if  tL  is  replaced  by  t^  (time  of  arrival)  and 
r(L)  is  replaced  by  r(A)  (position  at  tA). 
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Appendix  D 

Assumed  System  Errors 

This  appendix  contains  the  assumed  r.m.s.  injection  errors, 

the  initial  correlation  matrix  (EQ)  based  on  these  errors, 
and  the  assumed  variance  of  our  optical  measuring  device 
(Ref  2:28-29) . 

Initial  Correlation  Matrix  ( EQ) 

The  correlation  matrix  of  the  injection  errors  is  based  on 
the  assumed  r.m.s.  injection  errors  contained  in  Table  D-I  . 


Table  D-I 

Assumed  R.M.S.  Injection  Errors 


Altitude 

Track 

Range 

Position 

10,000  ft 

15.000  ft 

5,000  ft 

Velocity 

15  ft/sec 

6  ft/sec 

4  ft/ sec 

The  correlation  matrix  below  was  obtained  by  a  transformation 
from  the  altitude,  track,  and  range  coordinate  system  of  Table 
D-I  to  the  (a,b,c)  coordinate  system  defined  in  Fig.  C-l.  The 
basic  units  of  the  matrix  as  shewn  in  «vq.  B-l  are  lunar  units 
X  10~10  and  lunar  units  per  day  X  10"^. 
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1.087 

-1.019 

-Q.413 

0 

0 

0 

-1.019 

0.937 

0.388 

0 

0 

0 

-0.413 

0.388 

0.137 

0 

c 

0 

Eo  = 

0 

0 

0 

0.357 

-2.353 

-0.647 

0 

0 

0 

-2.553 

11.692 

2.963 

0 

0 

0 

-0,647 

2.963 

o.75i 

( D-l) 

The  initial  uncertainties  in  position  and  Telocity  are  3*343  ailea  and 
11*374  miles  per  hour. 


Variance(  (f  n2)  of  the  Optical  Sextant 

The  standard  deviation  of  the  optical  measuring  device  (Bextant)  was 
assumed  to  be  approximately  0.G5  milliradiana  (10  aeconda  of  arc).  However, 
since  all  the  angular  measurements  in  this  study  are  made  relatively  close 
to  the  earth  and  moon,  an  additional  factor  was  introduced  in  same  of  the 
computer  runs  to  accomodate  the  decrease  in  accuracy  as  the  spacecraft 
gets  close  to  a  planet.  Based  on  these  considerations  the  sextant  when 
measuring  the  angle  between  the  earth  and  moon  was  assumed  to  have  a 
random  error  whose  variance 

2  1  2 

(0.00003)^  ♦  (rftn))  rad,  if  r(tn)<(  Z(tn) 

,  (D-2) 

(0.00003 )2  +  (zf^))^  rad,  if  r(tn)>  2(1^) 

where  r(tn)  and  Z( tn)  are  the  distances  in  miles  fro*  the  spacecraft  to 

the  earth  and  moon,  respectively.  For  measurements  involving  only  one 

plai  ^  t ,  the  appropriate  variance  was  used ;  for  example ,  the  measurement 

of  the  angle  between  the  earth  and  a  star  would  use  only  the  first  part 

of  Eq.  D-2. 
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